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Abstract 

We consider a facility location problem, where the objective is to "dis- 
perse" a number of facilities, i.e., select a given number k of locations from 
a discrete set of n candidates, such that the average distance between 
selected locations is maximized. In particular, we present algorithmic re- 
sults for the case where vertices are represented by points in d-dimensional 
space, and edge weights correspond to rectilinear distances. Problems of 
this type have been considered before, with the best result being an ap- 
proximation algorithm with performance ratio 2. For the case where k is 
fixed, we establish a linear-time algorithm that finds an optimal solution. 
For the case where k is part of the input, we present a polynomial-time 
approximation scheme. 

1 Introduction 

A common problem in the area of facility location is the selection of a given 
number of k locations from a set P of n feasible positions, such that the se- 
lected set has optimal distance properties. Natural objective functions are the 
maximization of the minimum distance, or of the average distance between se- 
lected points; dispersion problems of this type come into play whenever we want 
to minimize interference between the corresponding facilities. Examples include 
oil storage tanks, ammunition dumps, nuclear power plants, hazardous waste 
sites, and fast-food outlets (see In the latter paper, the problem of 

maximizing the average distance is called the Remote Clique problem. 

Formally, problems of this type can be described as follows: given a graph 
G = (V 7 , E) with n vertices, and non-negative edge weights w VliV2 — d(vx,V2) for 
(vx, v%) S E. Given k G {2, . . . , n}, find a subset S C V with 15*1 = k, such that 
W (S) '■— J2(v v )eE(S) d(vi,Vj) is maximized. (Here, E(S) denotes the edge set 
of the subgraph of G induced by the vertex set S.) 

*A preliminary version of this paper appears in the proceedings of APPROX 2000, pp. 
132-143. 
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From a graph theoretic point of view, this problem has been called a heaviest 
subgraph problem. Being a weighted version of a generalization of the problem of 
deciding the existence of a k- clique, i.e., a complete subgraph with k vertices, the 
problem is strongly NP-hard [T2|. It should be noted that Hastad showed 
that the problem Clique of maximizing the cardinality of a set of vertices with 
a maximum possible number of edges is in general hard to approximate within 
n 1 ^ 6 . For the heaviest subgraph problem, we want to maximize the number of 
edges for a set of vertices of given cardinality, so Hastad's result does not imply 
an immediate performance bound. 

Related Work 

Over recent years, there have been a number of approximation algorithms 
for various subproblems of this type. Feige and Seltser [7] have studied the 
graph problem (i.e., edge weights are or 1) and showed how to find in time 

n O((l+logf)/ £ ) a fc _ get g c y with w ^ > ^ _ £ ) / k\ providcd that a 

fc-clique exists. They also gave evidence that for k ~ n 1 / 3 , semidefinite pro- 
gramming fails to distinguish between graphs that have a /c-clique, and graphs 
with densest fc-subgraphs having average degree less than logrt. 

Kortsarz and Peleg ^2] describe a polynomial algorithm with performance 
guarantee O(n - 3885 ) for the general case where edge weights do not have to obey 
the triangle inequality. A newer algorithm by Feige, Kortsarz, and Peleg [5] gives 
an approximation ratio of 0(n 1/ ' 3 logn). For the case where k = £l(n), Asahiro, 
Iwama, Tamaki, and Tokuyama |3] give a greedy constant factor approximation, 
while Srivastav and Wolf |15| use semidefinite programming for improved per- 
formance bounds. For the case of dense graphs (i.e., \E\ — fi(n 2 )) and k = f2(n), 
Arora, Karger, and Karpinski JI] give a polynomial time approximation scheme. 
On the other hand, Asahiro, Hassin, and Iwama [2| show that deciding the ex- 
istence of a "slightly dense" subgraph, i.e., an induced subgraph on k vertices 
that has at least fl(k 1+6 ) edges, is NP-complete. They also showed it is NP- 
complete to decide whether a graph with e edges has an induced subgraph on k 
vertices that has ^-(1 + 0(w e_1 )) edges; the latter is only slightly larger than 

— vk-k ) ' wn i cn ls the the average number of edges in a subgraph with k 
vertices. 

For the case where edge weights fulfill the triangle inequality, Ravi, Rosen- 
krantz, and Tayi |14| give a heuristic with time complexity 0(n 2 ) and prove 
that it guarantees a performance bound of 4. (See Tamir |17j with reference to 
this paper.) Hassin, Rubinstein, and Tamir |10| give a different heuristic with 
time complexity 0(n 2 + k 2 log fc) with performance bound 2. On a related note, 
see Chandra and Halldorsson [S], who study a number of different remoteness 
measures for the subset k, including total edge weight w(S). If the graph from 
which a subset of size k is to be selected is a tree, Tamir shows that an 
optimal weight subset can be determined in O(nk) time. 
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In many important cases there is even more known about the set of edge 
weights than just the validity of triangle inequality This is the case when the 
vertex set V corresponds to a point set P in geometric space, and distances 
between vertices are induced by geometric distances between points. Given the 
practical motivation for considering the problem, it is quite natural to consider 
geometric instances of this type. In fact, it was shown by Ravi, Rosenkrantz, 
and Tayi in ^3] that for the case of Euclidean distances in two-dimensional 
space, it is possible to achieve performance bounds that are arbitrarily close 
to ir/2 ps 1.57 For other metrics, however, the best performance guarantee is 
the factor 2 by (TJJ- Despite of these approximation results, it should be noted 
that the complexity status of the problem is still open, i.e., it is known known 
whether the problem is NP-hard. 

An important application of our problem is data sampling and clustering, 
where points are to be selected from a large more-dimensional set. Different 
metric dimensions of a data point describe different metric properties of a cor- 
responding item. Since these properties are not geometrically related, distances 
are typically not evaluated by Euclidean distances. Instead, some weighted L\ 
metric is used. (See Erkut |H|.) For data sampling, a set of points is to be 
selected that has high average distance. For clustering, a given set of points 
is to be subdivided into k clusters, such that points from the same cluster are 
close together, while points from different clusters are far apart. If we do the 
clustering by choosing k center points, and assigning any point to its nearest 
cluster center, we have to consider the same problem of finding a set of center 
points with large average distance, which is equivalent to finding a fc-clique with 
maximum total edge weight. 

For results on maximizing the minimum L\ distance within a selected set of 
n points see Baur and Fekete 0], who showed that finding such a set within a 
given polygon cannot be approximated arbitrarily well, unless P=NP. Finally, 
Gritzmann, Klee, and Larmann [9] have studied a somewhat related geometric 
selection problem: Given a set v of m points in n-dimensional space, choose a 
subset of n + 1 points, such that the total volume of the resulting simplex is 
maximum. They showed that this problem is NP-hard when n is part of the 
input. 

Main Results 

In this paper, we consider point sets P in <i-dimensional space, where d is some 
constant. For the most part, distances are measured according to the rectilinear 
"Manhattan" norm L\. 

Our results include the following: 

• A linear time (O(n)) algorithm to solve the problem to optimality in case 
where k is some fixed constant. This is in contrast to the case of Euclidean 
distances, where there is a well-known lower bound of Q(n\ogn) in the 
computation tree model for determining the diameter of a planar point 
set, i.e., the special case d = 2 and k = 2 (see ^3]). 
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• A polynomial time approximation scheme for the case where k is not fixed. 
This method can be applied for arbitrary fixed dimension d. For the case 
of Euclidean distances in two-dimensional space, it implies a performance 
bound of y/2 + e, for any given e > 0. 

2 Preliminaries 

For the most part of this paper, all points are assumed to be points in the plane. 
Spaces of arbitrary fixed dimension will be discussed in the end. Distances are 
measured using the L\ norm, unless noted otherwise. The x- and y-coordinates 
of a point p are denoted by x p and y p . If p and q are two points in the plane, 
then the distance between p and q is d(p, q) = \x p — x q \ + \y p — y q \. We say that 
q is above p in direction of a vector c, if the inner products satisfy (q, c) > (p, c) 
We say that a point p is maximal in direction c with respect to a set of points 
P if it maximizes the inner product {(c,x) | x € P}. For example, if p is an 
element of a set of points P and p has a maximal y-coordinate, then p is maximal 
in direction (0,1) with respect to P, and a point p with minimal x-coordinate 
is maximal in direction (-1,0) with respect to P. If the set P is clear from the 
context, we simply state that p is maximal in direction c. 

The weight of a set of points P is the sum of the distances between all pairs 
of points in this set, and is denoted by w(P). Similarly, w(P,Q) denotes the 
total sum of distances between two sets P and Q. For L\ distances, w x (P) and 
w x (P, Q) denote the sum of x-distances within P, or between P and Q. 

3 Cliques of Fixed Size 

Let S = {so, si, • • • , Sfc-i} be a maximum weight subset of P, where k is a fixed 
integer greater than 1. We will label the x- and y-coordinates of a point seS 
by some (x a ,yb) with < a < k and < b < k such that x < x\ < . . . < Xk~i 
and yo <y\< ■ ■ ■ < Vk-i- (Note that in general, a ^ b for a point s = (x a ,yb)-) 
Then 

0<i<j<k 0<i<j<k 

Now we can use local optimality to reduce the family of subsets that we need 
to consider: 

Lemma 1 There is a maximum weight subset S' of P of cardinality k, such 
that each point in S' is maximal in direction (2i + 1 — k, 2 j + 1 — k) with respect 
to P\S' for some values of i and j with < i, j < k. 

Proof: Consider a maximum weight subset S C P of cardinality k. Let 
s = (xi,yi>) be a point in S, such that there are i points s/ = (xi,yi>) € S\ {s} 
with xi < Xi (i.e., to the left of s) and k — i — 1 points s; = (xi,yi>) e S'\{s} with 
xi > Xi (i.e., strictly to the right of s). Similarly let there be j points below s 
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and k — j — 1 points strictly above s. We claim that s is maximal in direction 
(2i + 1 - fc, 2j + 1 - fc) with respect to P - S. 

Consider replacing s by a point s' = s + (xh,Uh) hi P — S. Let 5 — (2i + 
1 - k)x h + (2j + 1 - k)y h . Let S" = 5 \ {s} U {s'}. Assume first that point 
s' = {x[ , y[, ) is such that x[ and y^, have the same rank in S' as x\ and yi' have 
in S, i.e., there are i points si = (xi,yi) G S' \ {s 1 } with xi < x\ and j points 
s; = (x; , y;>) G S" \ {s'} with yi> < y' v . Replacing s by s' changes the a;-distances 
to the points left of s by ixh, and the x-distances to the points right of s by 
(k — i — 1)(— Xh)- Similarly, the y-distances change by jyu and (k — j — l)yh- 
So w(S') = w(S) + (2i + 1 - k)x h + (2j + 1 - k)y h = w(S) + 5. Since w(S) is 
maximum, we derive that S < 0, i.e., no point in P — S is above any point in S 
in direction (2i + 1 — k, 2 j + 1 — k). 

If the x- and y-coordinates of s' do not have the same rank in S' as the x- and 
y-coordinates of s in S, then it is not hard to show that w(S') > w(S)+S, so 5 < 
0. Therefore in this case, s' is strictly below s in direction (2i+ 1 — fc, 2j + 1 — fc). 

We can also conclude that if s = (a;, , y^/ ) and s' = (x[ , y^ ) are at the same 
level in direction (2i + 1 — k,2j + 1 — k), i.e., if <5 = 0, then the x- and y- 
coordinates of s' do have the same rank in 5" as the x- and y-coordinates of s 
in 5 and w(S) = w{S'). □ 

Theorem 1 Given a constant value for k, a maximum weight subset S of a set 
of n points P , such that S has cardinality k, can be found in linear time. 

Proof: Consider all directions of the form (2i+l— k, 2j + l — k) with < i, j < k. 
For each direction (a, 6), find Sk(a,b), a set of fc points that are maximal in 
direction (a, b) with respect to P — Sk(a, b). Compute the set USk(a, b) and try 
all possible subsets of size k of this set until a subset of maximum weight is 
found. 

Correctness follows from the fact that Lemmani m P nes that 5 C U5fc(a, b). 
Since k is a constant, each set Sk(a,b) can be found in linear time. Since the 
cardinality of U5fc(a, b) is less than or equal to A; 3 , the result follows. From the 
discussion at the end of the proof of Lemma^Jwe can conclude that if the set of 
fc points maximal in a direction (a, b) is not unique, any set of k points maximal 
in this direction will work equally well. 

□ 

Note that in the above estimate, we did not try to squeeze the constants in 
the 0(n) running time. A closer look shows that for k — 2, not more than 2 
subsets of P need to be evaluated for possible optimality, for k = 3, 8 subsets 
are sufficient. 
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4 Cliques of Variable Size 



In this section we consider the scenario where k is not fixed, i.e., k is part of the 
input. We show that there is a polynomial time approximation scheme (PTAS), 
i.e., for any fixed positive e, there is a polynomial approximation algorithm that 
finds a solution that is within (1 + e) of the optimum. 

The basic idea is to use for each of the d coordinates a suitable subset of 
m e coordinate values that subdivide an optimal solution into subsets of equal 
cardinality. More precisely, we describe the case d — 2; we find (by enumeration) 
a subdivision of an optimal solution into m £ x m e rectangular cells CV, , each 
of which must contain a specific number kij of selected points. From each cell 
Cij, the points are selected in a way that guarantees that the total distance 
to all other cells except for the m e — 1 cells in the same "horizontal" strip or 
the m s — 1 cells in the same "vertical" strip is maximized. As it turns out, 
this can be done in a way that the total neglected distance within the strips is 
bounded by a fraction of (5m e — 9)/(2(m e — 1)(to e — 2)) of the weight of an 
optimal solution, yielding the desired approximation property. See Figure [T] for 
the overall picture. 
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Figure 1: Subdividing the plane into cells 



For ease of presentation we assume that A: is a multiple of m e and m e > 
2. Approximation algorithms for other values of k can be constructed in a 
similar fashion. Consider an optimal solution of k points, denoted by OPT. 
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Furthermore consider a division of the plane by a set of m E + 1 x-coordinatcs 
Co < ••• < £1 < £m s - Let Xi := {p = (x,y) eft 2 \ & <x < & +1 ,0 < i < m E } be 
the vertical strip between coordinates & and By enumeration of possible 

choices of £o, • • ■ , £ro e we ma Y assume that the £i have the property that, for 
an optimal solution, from each of the m E strips Xi precisely k/m E points of P 
are chosen. (A small perturbation does not change optimality or approximation 
properties of solutions. This shows that in case of several points sharing the 
same coordinates, ties may be broken arbitrarily; in that case, points on the 
boundary between two strips may be considered belonging to one or the other 
of those strips, whatever is convenient to reach the appropriate number of points 
in a strip.) 

In a similar manner, suppose we know m E + 1 y-coordinates r/o < r\\ < . . . < 
Vm? such that from each horizontal strip Yi :— {p — (x,y) £ 5ft 2 | rji < y < 
Vi+ii < i < m E } a subset of k/m E points are chosen for an optimal solution. 

Let Cij :— Xi DYj, and let kij be the number of points in OPT that are 
chosen from Cy. Since J2o<Km E % = J2o<j<m E % = k/m £l we may assume 
by enumeration over the 0(k™) possible partitions of k/m £ into m £ pieces that 
we know all the numbers fcy . 

Finally, define the vector V,-j := ({2i + 1 — m E )k/m e , (2j + 1 — m E )k/m E ). 
Now our approximation algorithm is as follows: from each cell CV,-, choose some 
k^ points that are maximal in direction Vy. (Overlap between the selections 
from different cells is avoided by proceeding in lexicographic order of cells, and 
choosing the fey points among the candidates that are still unselected.) Let 
HEUbe the point set selected in this way. 

It is clear that HEU can be computed in polynomial time. We will proceed 
by a series of lemmas to determine how well w(HEU) approximates w(OPT). 
In the following, we consider the distances involving points from a particular 
cell C^. Let HEUij be the set of kij points that are selected from Cy by the 
heuristic, and let OPTij be a set of fey points of an optimal solution that are 
attributed to Cy. Let Sy = OPT^ n HEUij. Furthermore we define S'y = 
HEUij \ OPT^, and = OPT id \ HEU tj . Let HEU lt , OPTi., HEU. 3 and 
OPT m j be the set of k/m £ points selected from Xi and Yj by the heuristic and 
an optimal algorithm respectively. Finally HEU, :— HEU\ HEU,, HEU,j := 
HEU\ HEU.j, ~OPTi. := OPT\ OPTi. and OPT.j := OPT\ OPT.j. 

Lemma 2 



w x (HEUij,HEUi.) + w y {HEUij,HEU.j) 
> w^OPTij, OPTi.) +w y (OPTij, OPT.j). 

Proof: Consider a point p £ Sij. Thus, there is a point p' S S%j that was 
chosen by the heuristic instead of p. Now we can argue like in Lemma ^ Let 
h = (h x , h y ) = p'—p. When replacing p in OPT 'by p' , we increase the x-distance 
to the ik/m e points "left" of Cy by h x , while decreasing the x-distance to 
(m e — i—l)k/m s points "right" ofCy by h x . In the balance, this yields a change 
of ((2i + 1 — m £ )k/m e )h x . Similarly, we get a change of ((2j + 1 — m E )k/m e )h y 
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for the y-coordinates. By definition, we have assumed that the inner product 
(h, Vy) > 0, so the overall change of distances is nonnegative. 

Performing these replacements for all points in OPT\ HEU, we can trans- 
form OPT to HEU, while increasing the sum of distances w x (OPTij, OPTt,) + 
w y (OPT lj ,'OPT, j ) to the sum w x {HEU i] ,HEU l .) + w y (HEUij,HEU 9j ). 

□ 



In the following three lemmas we show that the total difference between the 
weight of an optimal solution w(OPT) and the total value of all the right hand 
sides (when summed over i) of the inequality in Lemma |3 is a small fraction of 
w(OPT). 

Lemma 3 

w x {OPT) 



w x {0PT i9 ) < 



2(m e - 2) 

Proof: Let Si — — Since i(m e —i—l)> m e — 2 for < i < m £ — 1, we 
have for < i < m e — 1 

w x (OPTi.) < — -^5i < 6 t — — . 

Zm £ m e m e Z(m £ — I) 

Since OPT has ik/m £ and (m e — i — l)k/m e points to the left of & and right of 
respectively, we have 



\ ik {m F — i — l)k „ 

■j x {OPT) > V — L_E >s, 

* — * m F r " - 



0<i<m E -l 
SO 

0<i<m E -l v ' 

□ 

Lemma 4 For i = and i = m e — 1 we have 

w x (OPT) 



w x {OPTi.) < 



m £ - 1 



Proof: Without loss of generality assume i = 0. Let xq, Xx, — , £(fc/m e )-i be 
the x-coordinates of the points po,p\, . ■ . ,P(k/m e )-i m OPT 0m . So 

fc fc 
w x (OPT .) = ( i - x ) + ( 3)(x_i^_ 2 - Xx) + ... 

< (— - i)(a - xo) + (— - s)(6 - + . . . 

k . . .. k . . . 

< — (a - a* + — (a - xi) + ••• 



Since £1 — Xj < x — Xj where < j < k/m e and x is the ^-coordinate of any 
point in OPTq, and since there are (m s — l)k/m E points in OPT^ mi we have 

fr-*i< (m J!: 1)fc ^fe^°-) 

so 



0<i<, 



< — !— r V iw^Cpi, OPT .) 

TO F — 1 * — ' 



1 



m £ — 1 



w x (OPT „ OPT .) 



< —^—w x (OPT). 

m e -l 

□ 

This proves the main properties. Now we only have to combine the above 
estimates to get an overall performance bound: 

Lemma 5 

£ w x {OPT l .,~OPT l .) + £ w v {OPT. h ~OPT. j ) 

0<i<m e 0<j<m e 

> (i _ — 5m -9 — -wopm 

" V 2(m E -l)(m E -2) ; V " 

Proof: From Lemmas (applied for indices < i < m — 1) and 01 (applied 
twice, once for i — 0, and once for i = m — 1), we derive that 

£ -^-) * 2(^ 5m r ) ( TO !-2) w -( o ^ 

0<i<m e V e A e ' 



and similarly 

.,/, (HPT -\ < 

2(m e -l)(m £ -2) 



0<i<m e 



Since 



w(OPT) = w x (OPT) + w y (OPT) 

£ w x (OPT l .,~dPT t .)+ £ ^(OPT 1 ,.) 

0<i<m e 0<i<m e 

+ £ w y (OPT. J ,OPT. J )+ £ w v( 0PT v)> 

0<j<m £ 0<j<m £ 
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the result follows. □ 



Putting together Lemma [21 and the error estimate from Lemma the ap- 
proximation theorem can now be proven. 

Theorem 2 For any fixed m, HEU can be computed in polynomial time, and 

^ S (1 - 2(„„ 5 '"lV-!-2) ""'° PT )^ 
The running time is exponential in i. 

Proof: The claim about the running time is clear. (The only step that is 
exponential in - is the enumeration over all 0(k me ) possible partitions of k/m E 
into m E pieces.) Using Lemmas El and El we derive 

w(HEU)> w x (HEUi.,HEUi.) + ^ w y (HEU. j ,HEU. j ) 

0<i<m e 0<j<m e 

> Yl MOPT l .,OPT l .)+ w y (OPT. j ,OPT. j ) 

> (1 - -, 5me ~ 9 -)w(OPT). 

~ K 2(m £ -l)(m £ -2) J v ' 

□ 



5 Implications 

It is straightforward to modify our above arguments to point sets under L\ 
distances in an arbitrary c?-dimensional space, with fixed d. 

Theorem 3 Given a constant value for k and d, a maximum weight subset S 
of a set of n points in d-dimensional space, such that S has cardinality k, can 
be found in linear time. If d and e are constants, but k is not fixed, then there 
is a polynomial time algorithm that finds a subset whose weight is within (1 + e) 
of the optimum. 

For the case of fixed k, it is straightforward to generalize the argument from 
Section |3 to see that there are at most (2k) d interesting directions to consider. 
For k being part of the input, the approximation scheme can be generalized in a 
straightforward manner by using an m e -subdivision in each coordinate direction. 
Again the complexity ends up being exponential in ~, as well as in d. 

For the case of distances in the plane, the results for Li distances can be 
applied by a standard argument: A rotation by 7r/4 transforms L m distances 
into L\ distances and vice versa. Furthermore, we can use the approximation 
scheme from the previous section to get a v2(l + e) approximation factor for 
the case of Euclidean distances in two-dimensional space, for any e > 0: In 
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polynomial time, find a fe-set Si such that Li(S) is within (1 + e) of an optimal 
solution OPTi with respect to Li distances. Let OPTq, be an optimal solution 
with respect to L2 distances. Then 

i2(s) £ ^■ (s,s 7i(r^ i ' <opr ' , -7i(r^ il<OPT2) 

£ 7WTT) L ' (fiPn) - 

and the claim follows. Similarly, any norm has its characteristic approximation 
factor p with respect to Li or distances; this factor immediately yields a 
(p + e)-approximation for geometric dispersion. 



6 Conclusions 

We have presented algorithms for geometric instances of the maximum weighted 
fc-clique problem. Our results give a dramatic improvement over the previous 
best approximation factor of 2 that was presented in ^D] for the case of general 
metric spaces. This underlines the observation that geometry can help to get 
better algorithms for problems from combinatorial optimization. 

Furthermore, the algorithms in |10j give better performance for Euclidean 
metric than for Manhattan distances. We correct this anomaly by showing that 
among problems involving geometric distances, the rectilinear metric may allow 
better algorithms than the Euclidean metric. 

It remains an interesting open problem to show NP-hardness of a geometric 
version of the problem for spaces of fixed dimension. In particular, the case 
of Manhattan distances in the plane may actually turn out to be polynomially 
solvable. 
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